Stick-slip friction and nucleation dynamics of ultra-thin liquid films 
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We develop the theory for stick-slip motion in ultra-thin liquid films confined between two moving 
atomically-flat surfaces. Our model is based on hydrodynamic equation for the flow coupled to 
the dynamic order parameter field describing the "shear melting and freezing" of the confined fluid. 
This model successfully accounts for observed phenomenology of friction in ultra-thin films, including 
periodic and chaotic sequences of slips, transitions from stick-slip motion to steady sliding. 
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The nature of sliding friction is a fundamental phys- 
ical problem of prime practical importance While 
the possibility to create low-friction surfaces and lubri- 
cant fluids has been ubiquitous for almost all engineering 
applications, it has become crucial for design of mod- 
ern micro-miniature devices such as information storage 
and micro-electro mechanical systems, where low friction 
without stick-slip (or interrupted) motion is necessary. 

Studies of friction between atomically flat mica sur- 
faces separated by the ultra-thin layer of fluid lubri- 
cant have revealed a striking phenomenon S: in a cer- 
tain range of experimental parameters the fluid exhibited 
solid-like properties, in particular, a critical yield stress 
leading to stick-slips similar to that in solid-on-solid dry 
friction process Q. This behavior was attributed to the 
confinement-induced freezing of the lubricant and its re- 
curring melting due to the increasing shear stress: as the 
fluid thickness is reduced to several molecular layers, it 
freezes, but when the shear stress exceeds some critical 
value, it melts (so-called "shear melting" effect). This be- 
havior was confirmed by molecular dynamics simulations 
HP] that indicated ordering of the fluid due to confine- 
ment by the walls. 

A quest for the quantitative description of the stick-slip 
lubricant dynamics motivated several theoretical works 
. The important step has been made by Carlson and 
Batista Q who proposed phenomenological constitutive 
relation connecting the frictional forces to velocity and 
coordinates via the order parameter-like state variable 
reflecting the degree of melting. This model successfully 
described some of the observed phenomenology of the ex- 
periment H and gave a new insight into the dynamics. 
Yet many important questions including the very mech- 
anism of the onset of the stick-slip remain unresolved. 

In this Letter we develop a theory of a stick-slip motion 
in an ultra-thin confined liquids based on the equation 
for the flow coupled to the equation for the order pa- 
rameter (OP) for the melting transition in the presence 
of the shear stress. We propose that the shear melting 
is controlled by the stress tensor rather than the sliding 
velocity as assumed in M . Making use of the generalized 



Lindemann criterion, we combine shear and thermody- 
namic melting within a unified description. Using this 
approach we describe the onset of the stick-slip motion as 
the function of the film thickness and determine the dy- 
namic phase diagram. We demonstrate that random nu- 
cleation of droplets of the fluid phase during the motion 
leads to irregular temporal distribution of slip events. 

Model. The flow of liquid lubricant has to satisfy the 
momentum conservation law: 
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where Vi is a component of the fluid velocity, cr^ is the 
stress tensor, D/Dt = d t + vV is the material derivative, 
and po is the density of fluid. Assuming incompressibility 
we set po = 1 and divv = 0. We further assume that the 
lubricant satisfies Maxwell-type stress-stain relation: 



dt<Jij + i/n, , = pi , 
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where t/y = dv.i/dxj +dvj/dxi is the shear strain rate, p, 
is the shear modulus, and rj is the shear stress relaxation 
rate. Thus, the stress-strain relation includes both vis- 
cous flow and elastic restoring forces. The conventional 
shear viscosity is defined as v = p/r/. To describe the 
dynamic phase transition between solid and fluid states 
we take into account that the stress relaxation rate r\ 
is itself a function of the physical state of the material 
quantified near the melting transition by the OP p which 
is defined in such a way that p = 1 corresponds to the 
solid state and p = to the liquid state. Physical inter- 
pretation of the OP for various systems can be different, 
but for crystalline solids p can be related to the disloca- 
tion density. We restrict ourselves to the simplest depen- 
dence of the stress relaxation rate on p: r\ = r]o(l — p), 
r/o = const. This choice assures that the Eq. (g) gives 
the standard Hook's law for the pure solid (p = 1) and 
the standard viscous stress-strain relation for the New- 
tonian fluid with p = 0. We assume that the OP obeys 
the standard Ginzburg-Landau equation 



tq dtp 



l 2 \7 2 p- 



p(l- P )(S-p) 



(3) 



1 



where tq and I are the characteristic time and length 
correspondingly; I should be of the order of the lattice 
constant a, and the characteristic time can be expressed 
through the sound velocity c s , tq ~ l/c s . Here 6 is 
the control parameter proportional to the temperature 
T . Since the melting transition in the lubricating layer 
occurs under the out-of-equilibrium conditions, it is char- 
acterized by two critical temperatures, T\, corresponding 
to an absolute instability of the overcooled liquid and 
T 2 , the stability limit of the solid phase The 
conventional thermodynamic melting temperature, T m , 
is confined between these limits: T\ < T m < T 2 . The 
parameter S is naturally expressed in the form 



6=(T-T 1 )/(T 2 -T 1 ) 
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Now we have to relate the solid instability temperature 
T 2 to the stress generated in the process of motion. To 
this end we notice first that in the absence of dynamic 
shear deformations the temperature T% can be estimated 
as T 2 ° = c^/ia 3 , where a is the lattice constant for the 
solid crystalline state, \i is the shear modulus (we write 
everything in a scalar form for simplicity), and c 2 L is the 
numerical factor ( "Lindemann number" ) . Assuming in- 
dependence of thermal fluctuations and shear, one can 
present a rms displacement field in a form: 
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where a = a xy is the shear stress, and the first term 
in the rhs stands for the thermal average displacement, 
while the second term expresses the shear-induced dis- 
placement field. In writing Eq. (^) we hypothesize that 
the solid phase instability can stem not only from the 
thermal fluctuations, but also from the shear, generated 
by the mutual motion of the solid surfaces confining the 
lubricant. At zero physical temperature the instability 
can be caused by this shear only, this concept of shear- 
induced melting generalizes the hypothesis of the dynamic 
disorder-driven melting introduced in earlier work [ fL2| . 
At the temperature T 2 the relation u 2 = c 2 a 2 holds, and 
we immediately obtain from Eq.([5|) 
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Substituting Eq. 
parameter 6 Eq. 



s) into the expression for the control 
) one derives in the first order 



6 = 6 Q + a 2 /a 2 



(7) 



w here 6 = (T - T X )/(T% - T x ) and a a = 
■J fi(T 2 — Ti)/aSo has a meaning of the yield shear stress. 

In the thin layer approximation we assume that the 
thickness of the lubricant layer h is small and neglect 
the dependence of the shear stress a on the transverse 
coordinate z. We can also neglect the dependence of a 
on the longitudinal coordinate x if the sample size L is 



not very large. Indeed, if the "acoustic shear time" r a = 
L/c s is much smaller then any characteristic time scale 
of the problem (e.g stick-slip time), the spatial variation 
of the stress is negligible, and the shear stress becomes 
the function of time only [lij . In this approximation after 
the integration over the area of the sample, Eq. (0) gives: 



at + ljiJ J ( l ~ P) dxdz = ^ 1 
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where V is the relative velocity of the upper plate with 
respect to the bottom. 

Now we can further simplify the OP dynamics. Since 
the walls favor the formation of the solid, the boundary 
conditions for OP read: p(0) = p(h) = 1, and the bulk 
variation of the OP are small as compared to 1. Let us 
seek the solution in a form 



p(x, z, t) = 1 — A(x, t) sinin/hz) 
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where A -C 1 is the slowly varying amplitude. Substitut- 
ing Eq. (||) into Eq. (||) and making use of the standard 
orthogonality procedure leads to 



A t = A x 
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where S — 5q + a 2 , and variable were rescaled as x/l 
x,t/r — > t,a xy /cro — > a. Accordingly, Eq. (|§|) yields: 



a t H -a / A(x)dx = —, 

ttL J n 
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where vq = fj,V/cro is the normalized pulling velocity. 
Note that in the most of the dynamic friction experi- 
ments the pulling velocity V does not coincide with the 
upper plate velocity V. In fact, the relation between the 
position of the upper plate x, friction force F, and the 
position of the spring (neglecting the mass of the spring 
and upper plate) is as follows: F — k(Vt — X), where X 
is the horizontal position of the upper plate and k is the 
stiffness of the spring. Since F ~ aL, and V = dX/dt, 
one immediately sees that exclusion of the plate velocity 
V results in the equation that is qualitatively the same 
as Eq. (^) but has renormalized vq and r/Q. Thus to 
simplify further discussion we will consider V = V . 

Stick-slip motion. First we discuss the spatially- 
uniform motion. In this case Eqs.([l0|),(pT|) become a 
pair of coupled ordinary differential equations (ODE). 
Above the melting temperature 6a > 1, the solid phase 
of the lubricant is formed due to proximity-to-the walls 
effect. In the experimentally relevant limit 7]q,Vq <C 1 
the abovementioned ODEs can be investigated analyt- 
ically by the bifurcation analysis and the multi-scale 
technique with a being a slow and A being a fast vari- 
able. Stick-slips are described by the limit cycle on the 
(7 — A plane. The transitions between different regimes 
are determined by the intersection of the manifold of 
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the slow motion -^-aA — Vo/h and the fast motion 

(S - 1 - f£) + ^(2 - 5) A - \A 2 = 0. The transition 
from the stick-slips to sliding corresponds to the inter- 
section of the slow motion manifold with the minimum 
of the fast motion manifold a — f(A), i.e. da/dA = 0. 
The limit cycle vanishes smoothly at the transition point 
if there is only one intersection of two manifolds, other- 
wise the cycle disappears abruptly with hysteresis. 
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FIG. 1. Phase of lubrication regimes at So = 1.3. Solid 
line indicates the continuous transition from stick-slip to slid- 
ing, and the dashed line corresponds to the hysteretic abrupt 
transition. Insets A,B show normalized shear stress a vs 
time for r/o = 0.01 at two different values of h correspond- 
ing to continuous and discontinuous transitions, and several 
values of vo/rjo approaching the transition line (dots in the 
main plot): h = 4, v /r]o = 0.1,0.2,0.3,0.31,0.32 (A) and 
h = 4.85 wo/770 = 0.03, 0.05, 0.064, 0.0645, 0.065 (B). 

Figures illustrate the transition from continuous 
sliding to the stick-slip motion. As one can see from 
Fig. [l], |^, stick-slips are possible only in relatively thin 
layers, while in the thick layers sliding is steady since 
the lubricant in the bulk is always in a fluid state. The 
critical thickness h c is determined by the stability con- 
dition da/dA = 0: h c = tt/^/Sq - 1 + 64(2 - S ) 2 /27tt 2 . 
We find that for h close to h c the transition from the 
stick-slip to sliding is always abrupt and has a hysteretic 
character (see Fig. || and Fig. |l|, Inset B). For the pa- 
rameters chosen, the "friction law" a vs vq has a mini- 
mum, (Fig. ^|, inset, curve 2) as common for a typical 
dry friction behavior Q . For smaller So the transition is 
continuous (inset A in Fig. [l|, curve 3 in Fig. 0, inset). 
For h < ho — n the dry friction without stick slip occurs 
because the viscous friction force become larger then the 
dry friction one, curve 1 in Inset | fL3f . 

Nucleation. The above analysis presumed that the 
stick-slip happens simultaneously in the entire space, 
which is certainly not the case for large samples. It is 
natural to expect that the stick-slip occurs via series of 



nucleation events when droplets of the liquid emerge in 
the solid phase and then expand and merge through- 
out the system. It is interesting to connect the nucle- 
ation dynamics in the systems with the lubricated fric- 
tion with the creation of topological defects during rapid 
quench ( "cosmological scenario") [l~5|-p^7|. Since in the 
stick phase the value of A is close to zero, the solid 
phase can be significantly "overheated" by the shear. For 
A -C 1, integration of Eq. ( |l(i| ) yields 

A & A exp (J (-e + a 2 )dt^j , (12) 

where e = l+fr-^o > 0. From Eq. (|l|), a w v t/h+a , 
and Aq,(Jo & r e the values of amplitude and stress in the 
beginning of stick phase. We restrict ourselves to the case 
h — > h c , where ag — > and A — 0(1). While a is small, 
the amplitude A decays exponentially to very small val- 
ues. It reaches a minimum A m i n — exp(— 2e 3 ^ 2 h/3vo) 
at t m in = € 1 ^ 2 h/vQ. Then it starts slowly to grow and 
reaches the value O(l) (slip event) at t = t m = y3t m i n . 
At t min < t < t m the growthrate — e + cr 2 > 0, and there- 
fore the lubricant is in an unstable ("overheated") state. 
The overheated solid is very sensitive to fluctuations, e.g. 
to thermal noise. Small nuclei of liquid can appear and 
expand within the solid, resulting to an accelerated slip 
event. Since at the low noise level nucleation events have 
the probabilistic character, one can expect the spatio- 
temporal randomness of the slip events. At the higher 
level of the noise the slips become more regular because 
the number of nucleation sites increases and the overall 
effect of the noise is averaged out. 
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FIG. 2. Temperature (<5o)-thickness (h) diagram. Above 
solid line (h c ) sliding occurs for arbitrary small velocity vq. 
Left of h c stick-slips exist for no < v c , while the transition 
to sliding is abrupt with hysteresis between solid and dashed 
lines and smooth otherwise. Below dotted line (ho) one has 
dry friction without stick-slips. Inset: shear stress a (or fric- 
tion force) vs sliding velocity vo in three different regions. 

We studied Eqs. ( ^0| , pTl ) numerically in a fairly large 
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domain, see Figs. ||,|]. Since during the slip phase the 
shear stress rapidly drops, the domains do not neces- 
sary propagate through the entire system and in large 
systems one may observe "partial slips." The random 
character of the nucleation process manifests itself in the 
non-periodicity of slips and local amplitude A. 




FIG. 3. Space-time plot of A for 8q = 1.1, h = 4, n = 0.01, 
v = 0.0002 in the system of length L = 1000. Black cor- 
respond to A — O(solid), white to A = 1 (liquid). Time 
progresses from top to bottom, total integration time 30000 
dimensionless units. Uncorrelated noise with zero average and 
amplitude 10 -16 is add on each time step and each grid point. 




FIG. 4. a and A at x = L/2 vs t for parameters of Fig. y. 

In conclusion, we have investigated the friction dynam- 
ics in ultra-thin liquid films. We have related the stick- 
slip behavior of thin liquid films to shear-induced melting 
and freezing processes. The developed approach allowed 
for the first time a quantitative description of dynamic 
nucleation effects leading to slip events. The proposed 



theory can be applied to a wide range of phenomena in- 
cluding friction in nanoscale devices, friction on ice, gran- 
ular materials [ [l8| as well as depinning transitions of flux- 
line lattices in type-II superconductors, charge density 
waves and other structures driven through disorder |l2[ . 
Our model offers a description for the long-standing prob- 
lem of the ultra-sound cmition during the friction dynam- 
ics. This research is supported by the Office of the Basic 
Energy Sciences at the US DOE, grants W-31-109-ENG- 
38, DE-FG03-95ER14516, and DE-FG03-96ER14592 
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